home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_tbmv.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  54.3 KB  |  1,820 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_tbmv (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int trans = 111;
  14.    int uplo = 121;
  15.    int diag = 131;
  16.    int N = 3;
  17.    int K = 1;
  18.    int lda = 3;
  19.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  20.    float X[] = { -0.089f, -0.688f, -0.203f };
  21.    int incX = -1;
  22.    float x_expected[] = { 0.017088f, 0.315595f, 0.243875f };
  23.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  24.    {
  25.      int i;
  26.      for (i = 0; i < 3; i++) {
  27.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 894)");
  28.      }
  29.    };
  30.   };
  31.  
  32.  
  33.   {
  34.    int order = 101;
  35.    int trans = 111;
  36.    int uplo = 121;
  37.    int diag = 132;
  38.    int N = 3;
  39.    int K = 1;
  40.    int lda = 3;
  41.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  42.    float X[] = { -0.089f, -0.688f, -0.203f };
  43.    int incX = -1;
  44.    float x_expected[] = { -0.089f, -0.721909f, 0.129992f };
  45.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  46.    {
  47.      int i;
  48.      for (i = 0; i < 3; i++) {
  49.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 895)");
  50.      }
  51.    };
  52.   };
  53.  
  54.  
  55.   {
  56.    int order = 101;
  57.    int trans = 111;
  58.    int uplo = 122;
  59.    int diag = 131;
  60.    int N = 3;
  61.    int K = 1;
  62.    int lda = 3;
  63.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  64.    float X[] = { -0.089f, -0.688f, -0.203f };
  65.    int incX = -1;
  66.    float x_expected[] = { 0.156927f, -0.159004f, 0.098252f };
  67.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  68.    {
  69.      int i;
  70.      for (i = 0; i < 3; i++) {
  71.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 896)");
  72.      }
  73.    };
  74.   };
  75.  
  76.  
  77.   {
  78.    int order = 101;
  79.    int trans = 111;
  80.    int uplo = 122;
  81.    int diag = 132;
  82.    int N = 3;
  83.    int K = 1;
  84.    int lda = 3;
  85.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  86.    float X[] = { -0.089f, -0.688f, -0.203f };
  87.    int incX = -1;
  88.    float x_expected[] = { 0.043096f, -0.584876f, -0.203f };
  89.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  90.    {
  91.      int i;
  92.      for (i = 0; i < 3; i++) {
  93.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 897)");
  94.      }
  95.    };
  96.   };
  97.  
  98.  
  99.   {
  100.    int order = 102;
  101.    int trans = 111;
  102.    int uplo = 121;
  103.    int diag = 131;
  104.    int N = 3;
  105.    int K = 1;
  106.    int lda = 3;
  107.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  108.    float X[] = { -0.089f, -0.688f, -0.203f };
  109.    int incX = -1;
  110.    float x_expected[] = { 0.024831f, -0.24504f, 0.447756f };
  111.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  112.    {
  113.      int i;
  114.      for (i = 0; i < 3; i++) {
  115.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 898)");
  116.      }
  117.    };
  118.   };
  119.  
  120.  
  121.   {
  122.    int order = 102;
  123.    int trans = 111;
  124.    int uplo = 121;
  125.    int diag = 132;
  126.    int N = 3;
  127.    int K = 1;
  128.    int lda = 3;
  129.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  130.    float X[] = { -0.089f, -0.688f, -0.203f };
  131.    int incX = -1;
  132.    float x_expected[] = { -0.089f, -0.670912f, 0.146504f };
  133.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  134.    {
  135.      int i;
  136.      for (i = 0; i < 3; i++) {
  137.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 899)");
  138.      }
  139.    };
  140.   };
  141.  
  142.  
  143.   {
  144.    int order = 102;
  145.    int trans = 111;
  146.    int uplo = 122;
  147.    int diag = 131;
  148.    int N = 3;
  149.    int K = 1;
  150.    int lda = 3;
  151.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  152.    float X[] = { -0.089f, -0.688f, -0.203f };
  153.    int incX = -1;
  154.    float x_expected[] = { -0.24504f, 0.447756f, -0.089117f };
  155.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  156.    {
  157.      int i;
  158.      for (i = 0; i < 3; i++) {
  159.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 900)");
  160.      }
  161.    };
  162.   };
  163.  
  164.  
  165.   {
  166.    int order = 102;
  167.    int trans = 111;
  168.    int uplo = 122;
  169.    int diag = 132;
  170.    int N = 3;
  171.    int K = 1;
  172.    int lda = 3;
  173.    float A[] = { 0.439f, -0.484f, -0.952f, -0.508f, 0.381f, -0.889f, -0.192f, -0.279f, -0.155f };
  174.    float X[] = { -0.089f, -0.688f, -0.203f };
  175.    int incX = -1;
  176.    float x_expected[] = { -0.351128f, -0.589748f, -0.203f };
  177.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  178.    {
  179.      int i;
  180.      for (i = 0; i < 3; i++) {
  181.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 901)");
  182.      }
  183.    };
  184.   };
  185.  
  186.  
  187.   {
  188.    int order = 101;
  189.    int trans = 112;
  190.    int uplo = 121;
  191.    int diag = 131;
  192.    int N = 3;
  193.    int K = 1;
  194.    int lda = 3;
  195.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  196.    float X[] = { 0.023f, -0.501f, -0.562f };
  197.    int incX = -1;
  198.    float x_expected[] = { 0.156047f, 0.189418f, -0.52828f };
  199.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  200.    {
  201.      int i;
  202.      for (i = 0; i < 3; i++) {
  203.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 902)");
  204.      }
  205.    };
  206.   };
  207.  
  208.  
  209.   {
  210.    int order = 101;
  211.    int trans = 112;
  212.    int uplo = 121;
  213.    int diag = 132;
  214.    int N = 3;
  215.    int K = 1;
  216.    int lda = 3;
  217.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  218.    float X[] = { 0.023f, -0.501f, -0.562f };
  219.    int incX = -1;
  220.    float x_expected[] = { 0.194342f, -0.449858f, -0.562f };
  221.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  222.    {
  223.      int i;
  224.      for (i = 0; i < 3; i++) {
  225.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 903)");
  226.      }
  227.    };
  228.   };
  229.  
  230.  
  231.   {
  232.    int order = 101;
  233.    int trans = 112;
  234.    int uplo = 122;
  235.    int diag = 131;
  236.    int N = 3;
  237.    int K = 1;
  238.    int lda = 3;
  239.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  240.    float X[] = { 0.023f, -0.501f, -0.562f };
  241.    int incX = -1;
  242.    float x_expected[] = { -0.0046f, 0.156047f, 0.189418f };
  243.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  244.    {
  245.      int i;
  246.      for (i = 0; i < 3; i++) {
  247.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 904)");
  248.      }
  249.    };
  250.   };
  251.  
  252.  
  253.   {
  254.    int order = 101;
  255.    int trans = 112;
  256.    int uplo = 122;
  257.    int diag = 132;
  258.    int N = 3;
  259.    int K = 1;
  260.    int lda = 3;
  261.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  262.    float X[] = { 0.023f, -0.501f, -0.562f };
  263.    int incX = -1;
  264.    float x_expected[] = { 0.023f, -0.516295f, -0.423724f };
  265.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  266.    {
  267.      int i;
  268.      for (i = 0; i < 3; i++) {
  269.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 905)");
  270.      }
  271.    };
  272.   };
  273.  
  274.  
  275.   {
  276.    int order = 102;
  277.    int trans = 112;
  278.    int uplo = 121;
  279.    int diag = 131;
  280.    int N = 3;
  281.    int K = 1;
  282.    int lda = 3;
  283.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  284.    float X[] = { 0.023f, -0.501f, -0.562f };
  285.    int incX = -1;
  286.    float x_expected[] = { 0.328565f, 0.326454f, 0.051142f };
  287.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  288.    {
  289.      int i;
  290.      for (i = 0; i < 3; i++) {
  291.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 906)");
  292.      }
  293.    };
  294.   };
  295.  
  296.  
  297.   {
  298.    int order = 102;
  299.    int trans = 112;
  300.    int uplo = 121;
  301.    int diag = 132;
  302.    int N = 3;
  303.    int K = 1;
  304.    int lda = 3;
  305.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  306.    float X[] = { 0.023f, -0.501f, -0.562f };
  307.    int incX = -1;
  308.    float x_expected[] = { 0.356165f, -0.345888f, -0.562f };
  309.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  310.    {
  311.      int i;
  312.      for (i = 0; i < 3; i++) {
  313.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 907)");
  314.      }
  315.    };
  316.   };
  317.  
  318.  
  319.   {
  320.    int order = 102;
  321.    int trans = 112;
  322.    int uplo = 122;
  323.    int diag = 131;
  324.    int N = 3;
  325.    int K = 1;
  326.    int lda = 3;
  327.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  328.    float X[] = { 0.023f, -0.501f, -0.562f };
  329.    int incX = -1;
  330.    float x_expected[] = { -0.015295f, 0.13041f, -0.482689f };
  331.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  332.    {
  333.      int i;
  334.      for (i = 0; i < 3; i++) {
  335.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 908)");
  336.      }
  337.    };
  338.   };
  339.  
  340.  
  341.   {
  342.    int order = 102;
  343.    int trans = 112;
  344.    int uplo = 122;
  345.    int diag = 132;
  346.    int N = 3;
  347.    int K = 1;
  348.    int lda = 3;
  349.    float A[] = { 0.94f, -0.091f, 0.984f, -0.276f, -0.342f, -0.484f, -0.665f, -0.2f, 0.349f };
  350.    float X[] = { 0.023f, -0.501f, -0.562f };
  351.    int incX = -1;
  352.    float x_expected[] = { 0.023f, -0.508866f, -0.516409f };
  353.    cblas_stbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  354.    {
  355.      int i;
  356.      for (i = 0; i < 3; i++) {
  357.        gsl_test_rel(X[i], x_expected[i], flteps, "stbmv(case 909)");
  358.      }
  359.    };
  360.   };
  361.  
  362.  
  363.   {
  364.    int order = 101;
  365.    int trans = 111;
  366.    int uplo = 121;
  367.    int diag = 131;
  368.    int N = 3;
  369.    int K = 1;
  370.    int lda = 3;
  371.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  372.    double X[] = { -0.77, -0.818, 0.337 };
  373.    int incX = -1;
  374.    double x_expected[] = { 0.50204, 0.563918, -0.590448 };
  375.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  376.    {
  377.      int i;
  378.      for (i = 0; i < 3; i++) {
  379.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 910)");
  380.      }
  381.    };
  382.   };
  383.  
  384.  
  385.   {
  386.    int order = 101;
  387.    int trans = 111;
  388.    int uplo = 121;
  389.    int diag = 132;
  390.    int N = 3;
  391.    int K = 1;
  392.    int lda = 3;
  393.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  394.    double X[] = { -0.77, -0.818, 0.337 };
  395.    int incX = -1;
  396.    double x_expected[] = { -0.77, -0.95429, -0.44419 };
  397.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  398.    {
  399.      int i;
  400.      for (i = 0; i < 3; i++) {
  401.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 911)");
  402.      }
  403.    };
  404.   };
  405.  
  406.  
  407.   {
  408.    int order = 101;
  409.    int trans = 111;
  410.    int uplo = 122;
  411.    int diag = 131;
  412.    int N = 3;
  413.    int K = 1;
  414.    int lda = 3;
  415.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  416.    double X[] = { -0.77, -0.818, 0.337 };
  417.    int incX = -1;
  418.    double x_expected[] = { 1.214016, -0.433258, 0.321835 };
  419.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  420.    {
  421.      int i;
  422.      for (i = 0; i < 3; i++) {
  423.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 912)");
  424.      }
  425.    };
  426.   };
  427.  
  428.  
  429.   {
  430.    int order = 101;
  431.    int trans = 111;
  432.    int uplo = 122;
  433.    int diag = 132;
  434.    int N = 3;
  435.    int K = 1;
  436.    int lda = 3;
  437.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  438.    double X[] = { -0.77, -0.818, 0.337 };
  439.    int incX = -1;
  440.    double x_expected[] = { -0.236664, -1.106472, 0.337 };
  441.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  442.    {
  443.      int i;
  444.      for (i = 0; i < 3; i++) {
  445.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 913)");
  446.      }
  447.    };
  448.   };
  449.  
  450.  
  451.   {
  452.    int order = 102;
  453.    int trans = 111;
  454.    int uplo = 121;
  455.    int diag = 131;
  456.    int N = 3;
  457.    int K = 1;
  458.    int lda = 3;
  459.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  460.    double X[] = { -0.77, -0.818, 0.337 };
  461.    int incX = -1;
  462.    double x_expected[] = { 0.68068, 0.357254, 1.022043 };
  463.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  464.    {
  465.      int i;
  466.      for (i = 0; i < 3; i++) {
  467.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 914)");
  468.      }
  469.    };
  470.   };
  471.  
  472.  
  473.   {
  474.    int order = 102;
  475.    int trans = 111;
  476.    int uplo = 121;
  477.    int diag = 132;
  478.    int N = 3;
  479.    int K = 1;
  480.    int lda = 3;
  481.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  482.    double X[] = { -0.77, -0.818, 0.337 };
  483.    int incX = -1;
  484.    double x_expected[] = { -0.77, -0.31596, 1.037208 };
  485.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  486.    {
  487.      int i;
  488.      for (i = 0; i < 3; i++) {
  489.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 915)");
  490.      }
  491.    };
  492.   };
  493.  
  494.  
  495.   {
  496.    int order = 102;
  497.    int trans = 111;
  498.    int uplo = 122;
  499.    int diag = 131;
  500.    int N = 3;
  501.    int K = 1;
  502.    int lda = 3;
  503.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  504.    double X[] = { -0.77, -0.818, 0.337 };
  505.    int incX = -1;
  506.    double x_expected[] = { 0.357254, 1.022043, 0.190742 };
  507.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  508.    {
  509.      int i;
  510.      for (i = 0; i < 3; i++) {
  511.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 916)");
  512.      }
  513.    };
  514.   };
  515.  
  516.  
  517.   {
  518.    int order = 102;
  519.    int trans = 111;
  520.    int uplo = 122;
  521.    int diag = 132;
  522.    int N = 3;
  523.    int K = 1;
  524.    int lda = 3;
  525.    double A[] = { 0.566, 0.955, -0.086, -0.856, 0.177, 0.974, -0.652, -0.884, 0.77 };
  526.    double X[] = { -0.77, -0.818, 0.337 };
  527.    int incX = -1;
  528.    double x_expected[] = { -0.914786, -0.496165, 0.337 };
  529.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  530.    {
  531.      int i;
  532.      for (i = 0; i < 3; i++) {
  533.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 917)");
  534.      }
  535.    };
  536.   };
  537.  
  538.  
  539.   {
  540.    int order = 101;
  541.    int trans = 112;
  542.    int uplo = 121;
  543.    int diag = 131;
  544.    int N = 3;
  545.    int K = 1;
  546.    int lda = 3;
  547.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  548.    double X[] = { -0.851, 0.481, 0.155 };
  549.    int incX = -1;
  550.    double x_expected[] = { 0.610833, -0.293243, 0.02914 };
  551.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  552.    {
  553.      int i;
  554.      for (i = 0; i < 3; i++) {
  555.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 918)");
  556.      }
  557.    };
  558.   };
  559.  
  560.  
  561.   {
  562.    int order = 101;
  563.    int trans = 112;
  564.    int uplo = 121;
  565.    int diag = 132;
  566.    int N = 3;
  567.    int K = 1;
  568.    int lda = 3;
  569.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  570.    double X[] = { -0.851, 0.481, 0.155 };
  571.    int incX = -1;
  572.    double x_expected[] = { -0.635031, 0.574, 0.155 };
  573.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  574.    {
  575.      int i;
  576.      for (i = 0; i < 3; i++) {
  577.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 919)");
  578.      }
  579.    };
  580.   };
  581.  
  582.  
  583.   {
  584.    int order = 101;
  585.    int trans = 112;
  586.    int uplo = 122;
  587.    int diag = 131;
  588.    int N = 3;
  589.    int K = 1;
  590.    int lda = 3;
  591.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  592.    double X[] = { -0.851, 0.481, 0.155 };
  593.    int incX = -1;
  594.    double x_expected[] = { 0.024679, 0.610833, -0.293243 };
  595.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  596.    {
  597.      int i;
  598.      for (i = 0; i < 3; i++) {
  599.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 920)");
  600.      }
  601.    };
  602.   };
  603.  
  604.  
  605.   {
  606.    int order = 101;
  607.    int trans = 112;
  608.    int uplo = 122;
  609.    int diag = 132;
  610.    int N = 3;
  611.    int K = 1;
  612.    int lda = 3;
  613.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  614.    double X[] = { -0.851, 0.481, 0.155 };
  615.    int incX = -1;
  616.    double x_expected[] = { -0.851, 0.875864, -0.231243 };
  617.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  618.    {
  619.      int i;
  620.      for (i = 0; i < 3; i++) {
  621.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 921)");
  622.      }
  623.    };
  624.   };
  625.  
  626.  
  627.   {
  628.    int order = 102;
  629.    int trans = 112;
  630.    int uplo = 121;
  631.    int diag = 131;
  632.    int N = 3;
  633.    int K = 1;
  634.    int lda = 3;
  635.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  636.    double X[] = { -0.851, 0.481, 0.155 };
  637.    int incX = -1;
  638.    double x_expected[] = { -0.198505, 0.091504, 0.093 };
  639.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  640.    {
  641.      int i;
  642.      for (i = 0; i < 3; i++) {
  643.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 922)");
  644.      }
  645.    };
  646.   };
  647.  
  648.  
  649.   {
  650.    int order = 102;
  651.    int trans = 112;
  652.    int uplo = 121;
  653.    int diag = 132;
  654.    int N = 3;
  655.    int K = 1;
  656.    int lda = 3;
  657.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  658.    double X[] = { -0.851, 0.481, 0.155 };
  659.    int incX = -1;
  660.    double x_expected[] = { -1.074184, 0.356535, 0.155 };
  661.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  662.    {
  663.      int i;
  664.      for (i = 0; i < 3; i++) {
  665.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 923)");
  666.      }
  667.    };
  668.   };
  669.  
  670.  
  671.   {
  672.    int order = 102;
  673.    int trans = 112;
  674.    int uplo = 122;
  675.    int diag = 131;
  676.    int N = 3;
  677.    int K = 1;
  678.    int lda = 3;
  679.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  680.    double X[] = { -0.851, 0.481, 0.155 };
  681.    int incX = -1;
  682.    double x_expected[] = { 0.394864, -0.768342, 0.31774 };
  683.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  684.    {
  685.      int i;
  686.      for (i = 0; i < 3; i++) {
  687.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 924)");
  688.      }
  689.    };
  690.   };
  691.  
  692.  
  693.   {
  694.    int order = 102;
  695.    int trans = 112;
  696.    int uplo = 122;
  697.    int diag = 132;
  698.    int N = 3;
  699.    int K = 1;
  700.    int lda = 3;
  701.    double A[] = { 0.188, 0.6, -0.743, -0.803, 0.449, -0.681, -0.464, -0.029, 0.553 };
  702.    double X[] = { -0.851, 0.481, 0.155 };
  703.    int incX = -1;
  704.    double x_expected[] = { -0.851, 0.098901, 0.4436 };
  705.    cblas_dtbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  706.    {
  707.      int i;
  708.      for (i = 0; i < 3; i++) {
  709.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtbmv(case 925)");
  710.      }
  711.    };
  712.   };
  713.  
  714.  
  715.   {
  716.    int order = 101;
  717.    int trans = 111;
  718.    int uplo = 121;
  719.    int diag = 131;
  720.    int N = 3;
  721.    int K = 1;
  722.    int lda = 3;
  723.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  724.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  725.    int incX = -1;
  726.    float x_expected[] = { -0.113114f, -0.051704f, -0.403567f, -0.288349f, -0.223936f, 0.841145f };
  727.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  728.    {
  729.      int i;
  730.      for (i = 0; i < 3; i++) {
  731.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 926) real");
  732.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 926) imag");
  733.      };
  734.    };
  735.   };
  736.  
  737.  
  738.   {
  739.    int order = 101;
  740.    int trans = 111;
  741.    int uplo = 121;
  742.    int diag = 132;
  743.    int N = 3;
  744.    int K = 1;
  745.    int lda = 3;
  746.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  747.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  748.    int incX = -1;
  749.    float x_expected[] = { -0.46f, 0.069f, -0.14027f, -0.23208f, -0.537722f, 0.841425f };
  750.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  751.    {
  752.      int i;
  753.      for (i = 0; i < 3; i++) {
  754.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 927) real");
  755.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 927) imag");
  756.      };
  757.    };
  758.   };
  759.  
  760.  
  761.   {
  762.    int order = 101;
  763.    int trans = 111;
  764.    int uplo = 122;
  765.    int diag = 131;
  766.    int N = 3;
  767.    int K = 1;
  768.    int lda = 3;
  769.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  770.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  771.    int incX = -1;
  772.    float x_expected[] = { -0.099689f, 0.487805f, 0.353793f, 0.325411f, -0.225658f, -0.776023f };
  773.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  774.    {
  775.      int i;
  776.      for (i = 0; i < 3; i++) {
  777.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 928) real");
  778.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 928) imag");
  779.      };
  780.    };
  781.   };
  782.  
  783.  
  784.   {
  785.    int order = 101;
  786.    int trans = 111;
  787.    int uplo = 122;
  788.    int diag = 132;
  789.    int N = 3;
  790.    int K = 1;
  791.    int lda = 3;
  792.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  793.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  794.    int incX = -1;
  795.    float x_expected[] = { -0.39057f, 0.113296f, 0.388863f, 0.131011f, -0.236f, 0.605f };
  796.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  797.    {
  798.      int i;
  799.      for (i = 0; i < 3; i++) {
  800.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 929) real");
  801.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 929) imag");
  802.      };
  803.    };
  804.   };
  805.  
  806.  
  807.   {
  808.    int order = 102;
  809.    int trans = 111;
  810.    int uplo = 121;
  811.    int diag = 131;
  812.    int N = 3;
  813.    int K = 1;
  814.    int lda = 3;
  815.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  816.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  817.    int incX = -1;
  818.    float x_expected[] = { -0.169119f, 0.443509f, 0.159816f, 0.139696f, -0.180955f, -0.835292f };
  819.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  820.    {
  821.      int i;
  822.      for (i = 0; i < 3; i++) {
  823.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 930) real");
  824.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 930) imag");
  825.      };
  826.    };
  827.   };
  828.  
  829.  
  830.   {
  831.    int order = 102;
  832.    int trans = 111;
  833.    int uplo = 121;
  834.    int diag = 132;
  835.    int N = 3;
  836.    int K = 1;
  837.    int lda = 3;
  838.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  839.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  840.    int incX = -1;
  841.    float x_expected[] = { -0.46f, 0.069f, 0.194886f, -0.054704f, -0.191297f, 0.545731f };
  842.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  843.    {
  844.      int i;
  845.      for (i = 0; i < 3; i++) {
  846.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 931) real");
  847.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 931) imag");
  848.      };
  849.    };
  850.   };
  851.  
  852.  
  853.   {
  854.    int order = 102;
  855.    int trans = 111;
  856.    int uplo = 122;
  857.    int diag = 131;
  858.    int N = 3;
  859.    int K = 1;
  860.    int lda = 3;
  861.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  862.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  863.    int incX = -1;
  864.    float x_expected[] = { 0.159816f, 0.139696f, -0.180955f, -0.835292f, 0.077786f, 0.60472f };
  865.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  866.    {
  867.      int i;
  868.      for (i = 0; i < 3; i++) {
  869.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 932) real");
  870.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 932) imag");
  871.      };
  872.    };
  873.   };
  874.  
  875.  
  876.   {
  877.    int order = 102;
  878.    int trans = 111;
  879.    int uplo = 122;
  880.    int diag = 132;
  881.    int N = 3;
  882.    int K = 1;
  883.    int lda = 3;
  884.    float A[] = { 0.824f, -0.45f, -0.987f, 0.758f, 0.42f, -0.357f, 0.147f, -0.191f, 0.88f, 0.63f, 0.155f, -0.573f, 0.224f, 0.146f, 0.501f, -0.889f, 0.456f, 0.796f };
  885.    float X[] = { -0.46f, 0.069f, 0.308f, -0.003f, -0.236f, 0.605f };
  886.    int incX = -1;
  887.    float x_expected[] = { -0.18707f, 0.2604f, 0.082342f, -0.779023f, -0.236f, 0.605f };
  888.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  889.    {
  890.      int i;
  891.      for (i = 0; i < 3; i++) {
  892.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 933) real");
  893.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 933) imag");
  894.      };
  895.    };
  896.   };
  897.  
  898.  
  899.   {
  900.    int order = 101;
  901.    int trans = 112;
  902.    int uplo = 121;
  903.    int diag = 131;
  904.    int N = 3;
  905.    int K = 1;
  906.    int lda = 3;
  907.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  908.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  909.    int incX = -1;
  910.    float x_expected[] = { 0.647885f, 0.621535f, -0.104407f, 0.05309f, 0.732704f, 0.055982f };
  911.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  912.    {
  913.      int i;
  914.      for (i = 0; i < 3; i++) {
  915.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 934) real");
  916.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 934) imag");
  917.      };
  918.    };
  919.   };
  920.  
  921.  
  922.   {
  923.    int order = 101;
  924.    int trans = 112;
  925.    int uplo = 121;
  926.    int diag = 132;
  927.    int N = 3;
  928.    int K = 1;
  929.    int lda = 3;
  930.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  931.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  932.    int incX = -1;
  933.    float x_expected[] = { 1.2955f, 0.190774f, -0.247934f, 0.982616f, -0.894f, -0.116f };
  934.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  935.    {
  936.      int i;
  937.      for (i = 0; i < 3; i++) {
  938.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 935) real");
  939.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 935) imag");
  940.      };
  941.    };
  942.   };
  943.  
  944.  
  945.   {
  946.    int order = 101;
  947.    int trans = 112;
  948.    int uplo = 122;
  949.    int diag = 131;
  950.    int N = 3;
  951.    int K = 1;
  952.    int lda = 3;
  953.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  954.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  955.    int incX = -1;
  956.    float x_expected[] = { 0.096482f, -0.071661f, 0.647885f, 0.621535f, -0.104407f, 0.05309f };
  957.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  958.    {
  959.      int i;
  960.      for (i = 0; i < 3; i++) {
  961.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 936) real");
  962.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 936) imag");
  963.      };
  964.    };
  965.   };
  966.  
  967.  
  968.   {
  969.    int order = 101;
  970.    int trans = 112;
  971.    int uplo = 122;
  972.    int diag = 132;
  973.    int N = 3;
  974.    int K = 1;
  975.    int lda = 3;
  976.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  977.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  978.    int incX = -1;
  979.    float x_expected[] = { 0.411f, -0.308f, -1.14861f, 0.933761f, -1.66247f, -0.234526f };
  980.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  981.    {
  982.      int i;
  983.      for (i = 0; i < 3; i++) {
  984.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 937) real");
  985.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 937) imag");
  986.      };
  987.    };
  988.   };
  989.  
  990.  
  991.   {
  992.    int order = 102;
  993.    int trans = 112;
  994.    int uplo = 121;
  995.    int diag = 131;
  996.    int N = 3;
  997.    int K = 1;
  998.    int lda = 3;
  999.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  1000.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  1001.    int incX = -1;
  1002.    float x_expected[] = { 0.632361f, -0.409373f, 0.578489f, 0.012724f, 0.664066f, 0.171616f };
  1003.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1004.    {
  1005.      int i;
  1006.      for (i = 0; i < 3; i++) {
  1007.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 938) real");
  1008.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 938) imag");
  1009.      };
  1010.    };
  1011.   };
  1012.  
  1013.  
  1014.   {
  1015.    int order = 102;
  1016.    int trans = 112;
  1017.    int uplo = 121;
  1018.    int diag = 132;
  1019.    int N = 3;
  1020.    int K = 1;
  1021.    int lda = 3;
  1022.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  1023.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  1024.    int incX = -1;
  1025.    float x_expected[] = { 0.946879f, -0.645712f, -1.21801f, 0.32495f, -0.894f, -0.116f };
  1026.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1027.    {
  1028.      int i;
  1029.      for (i = 0; i < 3; i++) {
  1030.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 939) real");
  1031.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 939) imag");
  1032.      };
  1033.    };
  1034.   };
  1035.  
  1036.  
  1037.   {
  1038.    int order = 102;
  1039.    int trans = 112;
  1040.    int uplo = 122;
  1041.    int diag = 131;
  1042.    int N = 3;
  1043.    int K = 1;
  1044.    int lda = 3;
  1045.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  1046.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  1047.    int incX = -1;
  1048.    float x_expected[] = { -0.236612f, 0.122761f, -1.12184f, -0.358823f, 1.4975f, -0.470595f };
  1049.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1050.    {
  1051.      int i;
  1052.      for (i = 0; i < 3; i++) {
  1053.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 940) real");
  1054.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 940) imag");
  1055.      };
  1056.    };
  1057.   };
  1058.  
  1059.  
  1060.   {
  1061.    int order = 102;
  1062.    int trans = 112;
  1063.    int uplo = 122;
  1064.    int diag = 132;
  1065.    int N = 3;
  1066.    int K = 1;
  1067.    int lda = 3;
  1068.    float A[] = { -0.814f, 0.043f, -0.755f, -0.094f, 0.876f, 0.257f, 0.406f, 0.491f, -0.27f, -0.787f, 0.545f, 0.732f, -0.512f, -0.085f, 0.234f, 0.001f, -0.225f, -0.002f };
  1069.    float X[] = { 0.411f, -0.308f, -0.912f, 0.811f, -0.894f, -0.116f };
  1070.    int incX = -1;
  1071.    float x_expected[] = { 0.411f, -0.308f, -1.26537f, 0.570703f, -0.129206f, -0.642577f };
  1072.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1073.    {
  1074.      int i;
  1075.      for (i = 0; i < 3; i++) {
  1076.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 941) real");
  1077.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 941) imag");
  1078.      };
  1079.    };
  1080.   };
  1081.  
  1082.  
  1083.   {
  1084.    int order = 101;
  1085.    int trans = 113;
  1086.    int uplo = 121;
  1087.    int diag = 131;
  1088.    int N = 3;
  1089.    int K = 1;
  1090.    int lda = 3;
  1091.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1092.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1093.    int incX = -1;
  1094.    float x_expected[] = { 0.413357f, 0.178267f, -0.114618f, -1.35595f, -0.513288f, 0.611332f };
  1095.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1096.    {
  1097.      int i;
  1098.      for (i = 0; i < 3; i++) {
  1099.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 942) real");
  1100.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 942) imag");
  1101.      };
  1102.    };
  1103.   };
  1104.  
  1105.  
  1106.   {
  1107.    int order = 101;
  1108.    int trans = 113;
  1109.    int uplo = 121;
  1110.    int diag = 132;
  1111.    int N = 3;
  1112.    int K = 1;
  1113.    int lda = 3;
  1114.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1115.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1116.    int incX = -1;
  1117.    float x_expected[] = { 0.368428f, 0.071217f, -0.954366f, -0.390486f, 0.694f, -0.954f };
  1118.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1119.    {
  1120.      int i;
  1121.      for (i = 0; i < 3; i++) {
  1122.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 943) real");
  1123.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 943) imag");
  1124.      };
  1125.    };
  1126.   };
  1127.  
  1128.  
  1129.   {
  1130.    int order = 101;
  1131.    int trans = 113;
  1132.    int uplo = 122;
  1133.    int diag = 131;
  1134.    int N = 3;
  1135.    int K = 1;
  1136.    int lda = 3;
  1137.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1138.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1139.    int incX = -1;
  1140.    float x_expected[] = { -0.084786f, -0.059464f, 0.413357f, 0.178267f, -0.114618f, -1.35595f };
  1141.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1142.    {
  1143.      int i;
  1144.      for (i = 0; i < 3; i++) {
  1145.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 944) real");
  1146.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 944) imag");
  1147.      };
  1148.    };
  1149.   };
  1150.  
  1151.  
  1152.   {
  1153.    int order = 101;
  1154.    int trans = 113;
  1155.    int uplo = 122;
  1156.    int diag = 132;
  1157.    int N = 3;
  1158.    int K = 1;
  1159.    int lda = 3;
  1160.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1161.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1162.    int incX = -1;
  1163.    float x_expected[] = { 0.065f, -0.082f, -0.636071f, 0.80005f, 0.787748f, -1.14446f };
  1164.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1165.    {
  1166.      int i;
  1167.      for (i = 0; i < 3; i++) {
  1168.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 945) real");
  1169.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 945) imag");
  1170.      };
  1171.    };
  1172.   };
  1173.  
  1174.  
  1175.   {
  1176.    int order = 102;
  1177.    int trans = 113;
  1178.    int uplo = 121;
  1179.    int diag = 131;
  1180.    int N = 3;
  1181.    int K = 1;
  1182.    int lda = 3;
  1183.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1184.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1185.    int incX = -1;
  1186.    float x_expected[] = { -1.18498f, -0.424201f, 0.230196f, 0.374209f, -0.208366f, -1.16549f };
  1187.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1188.    {
  1189.      int i;
  1190.      for (i = 0; i < 3; i++) {
  1191.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 946) real");
  1192.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 946) imag");
  1193.      };
  1194.    };
  1195.   };
  1196.  
  1197.  
  1198.   {
  1199.    int order = 102;
  1200.    int trans = 113;
  1201.    int uplo = 121;
  1202.    int diag = 132;
  1203.    int N = 3;
  1204.    int K = 1;
  1205.    int lda = 3;
  1206.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1207.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1208.    int incX = -1;
  1209.    float x_expected[] = { -1.03519f, -0.446737f, -0.819232f, 0.995992f, 0.694f, -0.954f };
  1210.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1211.    {
  1212.      int i;
  1213.      for (i = 0; i < 3; i++) {
  1214.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 947) real");
  1215.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 947) imag");
  1216.      };
  1217.    };
  1218.   };
  1219.  
  1220.  
  1221.   {
  1222.    int order = 102;
  1223.    int trans = 113;
  1224.    int uplo = 122;
  1225.    int diag = 131;
  1226.    int N = 3;
  1227.    int K = 1;
  1228.    int lda = 3;
  1229.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1230.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1231.    int incX = -1;
  1232.    float x_expected[] = { 0.109929f, 0.02505f, 0.062939f, -0.202464f, -0.470658f, 1.69006f };
  1233.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1234.    {
  1235.      int i;
  1236.      for (i = 0; i < 3; i++) {
  1237.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 948) real");
  1238.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 948) imag");
  1239.      };
  1240.    };
  1241.   };
  1242.  
  1243.  
  1244.   {
  1245.    int order = 102;
  1246.    int trans = 113;
  1247.    int uplo = 122;
  1248.    int diag = 132;
  1249.    int N = 3;
  1250.    int K = 1;
  1251.    int lda = 3;
  1252.    float A[] = { -0.675f, 0.047f, 0.695f, 0.724f, -0.438f, 0.991f, -0.188f, -0.06f, -0.093f, 0.302f, 0.842f, -0.753f, 0.465f, -0.972f, -0.058f, 0.988f, 0.093f, 0.164f };
  1253.    float X[] = { 0.065f, -0.082f, -0.746f, 0.775f, 0.694f, -0.954f };
  1254.    int incX = -1;
  1255.    float x_expected[] = { 0.065f, -0.082f, -0.776809f, 0.762996f, 0.73663f, 0.124729f };
  1256.    cblas_ctbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1257.    {
  1258.      int i;
  1259.      for (i = 0; i < 3; i++) {
  1260.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctbmv(case 949) real");
  1261.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctbmv(case 949) imag");
  1262.      };
  1263.    };
  1264.   };
  1265.  
  1266.  
  1267.   {
  1268.    int order = 101;
  1269.    int trans = 111;
  1270.    int uplo = 121;
  1271.    int diag = 131;
  1272.    int N = 3;
  1273.    int K = 1;
  1274.    int lda = 3;
  1275.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1276.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1277.    int incX = -1;
  1278.    double x_expected[] = { -0.010019, -0.1678, -0.042017, -1.112094, 0.010004, -0.480427 };
  1279.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1280.    {
  1281.      int i;
  1282.      for (i = 0; i < 3; i++) {
  1283.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 950) real");
  1284.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 950) imag");
  1285.      };
  1286.    };
  1287.   };
  1288.  
  1289.  
  1290.   {
  1291.    int order = 101;
  1292.    int trans = 111;
  1293.    int uplo = 121;
  1294.    int diag = 132;
  1295.    int N = 3;
  1296.    int K = 1;
  1297.    int lda = 3;
  1298.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1299.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1300.    int incX = -1;
  1301.    double x_expected[] = { 0.064, 0.169, -0.80842, -0.715637, -0.829924, -0.212971 };
  1302.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1303.    {
  1304.      int i;
  1305.      for (i = 0; i < 3; i++) {
  1306.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 951) real");
  1307.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 951) imag");
  1308.      };
  1309.    };
  1310.   };
  1311.  
  1312.  
  1313.   {
  1314.    int order = 101;
  1315.    int trans = 111;
  1316.    int uplo = 122;
  1317.    int diag = 131;
  1318.    int N = 3;
  1319.    int K = 1;
  1320.    int lda = 3;
  1321.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1322.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1323.    int incX = -1;
  1324.    double x_expected[] = { 0.634014, 0.796937, -0.585538, -0.895375, -0.125887, 0.010019 };
  1325.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1326.    {
  1327.      int i;
  1328.      for (i = 0; i < 3; i++) {
  1329.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 952) real");
  1330.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 952) imag");
  1331.      };
  1332.    };
  1333.   };
  1334.  
  1335.  
  1336.   {
  1337.    int order = 101;
  1338.    int trans = 111;
  1339.    int uplo = 122;
  1340.    int diag = 132;
  1341.    int N = 3;
  1342.    int K = 1;
  1343.    int lda = 3;
  1344.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1345.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1346.    int incX = -1;
  1347.    double x_expected[] = { 0.567497, 1.085122, -1.217792, -1.322566, -0.641, -0.103 };
  1348.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1349.    {
  1350.      int i;
  1351.      for (i = 0; i < 3; i++) {
  1352.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 953) real");
  1353.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 953) imag");
  1354.      };
  1355.    };
  1356.   };
  1357.  
  1358.  
  1359.   {
  1360.    int order = 102;
  1361.    int trans = 111;
  1362.    int uplo = 121;
  1363.    int diag = 131;
  1364.    int N = 3;
  1365.    int K = 1;
  1366.    int lda = 3;
  1367.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1368.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1369.    int incX = -1;
  1370.    double x_expected[] = { 0.130517, -0.119185, -0.187765, -0.519609, -0.169484, -1.165438 };
  1371.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1372.    {
  1373.      int i;
  1374.      for (i = 0; i < 3; i++) {
  1375.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 954) real");
  1376.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 954) imag");
  1377.      };
  1378.    };
  1379.   };
  1380.  
  1381.  
  1382.   {
  1383.    int order = 102;
  1384.    int trans = 111;
  1385.    int uplo = 121;
  1386.    int diag = 132;
  1387.    int N = 3;
  1388.    int K = 1;
  1389.    int lda = 3;
  1390.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1391.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1392.    int incX = -1;
  1393.    double x_expected[] = { 0.064, 0.169, -0.820019, -0.9468, -0.684597, -1.278457 };
  1394.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1395.    {
  1396.      int i;
  1397.      for (i = 0; i < 3; i++) {
  1398.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 955) real");
  1399.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 955) imag");
  1400.      };
  1401.    };
  1402.   };
  1403.  
  1404.  
  1405.   {
  1406.    int order = 102;
  1407.    int trans = 111;
  1408.    int uplo = 122;
  1409.    int diag = 131;
  1410.    int N = 3;
  1411.    int K = 1;
  1412.    int lda = 3;
  1413.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1414.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1415.    int incX = -1;
  1416.    double x_expected[] = { -0.187765, -0.519609, -0.169484, -1.165438, 0.198928, -0.370456 };
  1417.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1418.    {
  1419.      int i;
  1420.      for (i = 0; i < 3; i++) {
  1421.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 956) real");
  1422.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 956) imag");
  1423.      };
  1424.    };
  1425.   };
  1426.  
  1427.  
  1428.   {
  1429.    int order = 102;
  1430.    int trans = 111;
  1431.    int uplo = 122;
  1432.    int diag = 132;
  1433.    int N = 3;
  1434.    int K = 1;
  1435.    int lda = 3;
  1436.    double A[] = { -0.212, 0.612, 0.189, -0.046, -0.124, 0.82, 0.753, 0.727, 0.331, 0.116, 0.504, -0.673, -0.888, -0.277, -0.361, -0.909, 0.982, -0.124 };
  1437.    double X[] = { 0.064, 0.169, -0.81, -0.779, -0.641, -0.103 };
  1438.    int incX = -1;
  1439.    double x_expected[] = { -0.113746, -0.182809, -0.935887, -0.768981, -0.641, -0.103 };
  1440.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1441.    {
  1442.      int i;
  1443.      for (i = 0; i < 3; i++) {
  1444.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 957) real");
  1445.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 957) imag");
  1446.      };
  1447.    };
  1448.   };
  1449.  
  1450.  
  1451.   {
  1452.    int order = 101;
  1453.    int trans = 112;
  1454.    int uplo = 121;
  1455.    int diag = 131;
  1456.    int N = 3;
  1457.    int K = 1;
  1458.    int lda = 3;
  1459.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1460.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1461.    int incX = -1;
  1462.    double x_expected[] = { -0.436746, 0.963714, -1.087615, -0.018695, 0.30063, 0.12958 };
  1463.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1464.    {
  1465.      int i;
  1466.      for (i = 0; i < 3; i++) {
  1467.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 958) real");
  1468.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 958) imag");
  1469.      };
  1470.    };
  1471.   };
  1472.  
  1473.  
  1474.   {
  1475.    int order = 101;
  1476.    int trans = 112;
  1477.    int uplo = 121;
  1478.    int diag = 132;
  1479.    int N = 3;
  1480.    int K = 1;
  1481.    int lda = 3;
  1482.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1483.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1484.    int incX = -1;
  1485.    double x_expected[] = { 0.895682, 1.407174, 0.2408, -0.14282, -0.649, 0.188 };
  1486.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1487.    {
  1488.      int i;
  1489.      for (i = 0; i < 3; i++) {
  1490.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 959) real");
  1491.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 959) imag");
  1492.      };
  1493.    };
  1494.   };
  1495.  
  1496.  
  1497.   {
  1498.    int order = 101;
  1499.    int trans = 112;
  1500.    int uplo = 122;
  1501.    int diag = 131;
  1502.    int N = 3;
  1503.    int K = 1;
  1504.    int lda = 3;
  1505.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1506.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1507.    int incX = -1;
  1508.    double x_expected[] = { 0.785744, -0.3966, -0.436746, 0.963714, -1.087615, -0.018695 };
  1509.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1510.    {
  1511.      int i;
  1512.      for (i = 0; i < 3; i++) {
  1513.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 960) real");
  1514.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 960) imag");
  1515.      };
  1516.    };
  1517.   };
  1518.  
  1519.  
  1520.   {
  1521.    int order = 101;
  1522.    int trans = 112;
  1523.    int uplo = 122;
  1524.    int diag = 132;
  1525.    int N = 3;
  1526.    int K = 1;
  1527.    int lda = 3;
  1528.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1529.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1530.    int incX = -1;
  1531.    double x_expected[] = { 0.884, 0.636, 0.472572, 0.47454, -1.056415, 0.594125 };
  1532.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1533.    {
  1534.      int i;
  1535.      for (i = 0; i < 3; i++) {
  1536.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 961) real");
  1537.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 961) imag");
  1538.      };
  1539.    };
  1540.   };
  1541.  
  1542.  
  1543.   {
  1544.    int order = 102;
  1545.    int trans = 112;
  1546.    int uplo = 121;
  1547.    int diag = 131;
  1548.    int N = 3;
  1549.    int K = 1;
  1550.    int lda = 3;
  1551.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1552.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1553.    int incX = -1;
  1554.    double x_expected[] = { 0.464705, -0.108078, 0.094975, 0.376323, -0.6802, -0.42482 };
  1555.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1556.    {
  1557.      int i;
  1558.      for (i = 0; i < 3; i++) {
  1559.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 962) real");
  1560.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 962) imag");
  1561.      };
  1562.    };
  1563.   };
  1564.  
  1565.  
  1566.   {
  1567.    int order = 102;
  1568.    int trans = 112;
  1569.    int uplo = 121;
  1570.    int diag = 132;
  1571.    int N = 3;
  1572.    int K = 1;
  1573.    int lda = 3;
  1574.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1575.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1576.    int incX = -1;
  1577.    double x_expected[] = { 0.562961, 0.924522, 1.004293, -0.112851, -0.649, 0.188 };
  1578.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1579.    {
  1580.      int i;
  1581.      for (i = 0; i < 3; i++) {
  1582.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 963) real");
  1583.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 963) imag");
  1584.      };
  1585.    };
  1586.   };
  1587.  
  1588.  
  1589.   {
  1590.    int order = 102;
  1591.    int trans = 112;
  1592.    int uplo = 122;
  1593.    int diag = 131;
  1594.    int N = 3;
  1595.    int K = 1;
  1596.    int lda = 3;
  1597.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1598.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1599.    int incX = -1;
  1600.    double x_expected[] = { -0.448428, 0.19254, -0.674583, 1.236189, 0.780774, 1.167088 };
  1601.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1602.    {
  1603.      int i;
  1604.      for (i = 0; i < 3; i++) {
  1605.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 964) real");
  1606.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 964) imag");
  1607.      };
  1608.    };
  1609.   };
  1610.  
  1611.  
  1612.   {
  1613.    int order = 102;
  1614.    int trans = 112;
  1615.    int uplo = 122;
  1616.    int diag = 132;
  1617.    int N = 3;
  1618.    int K = 1;
  1619.    int lda = 3;
  1620.    double A[] = { -0.374, -0.308, 0.792, 0.884, -0.794, -0.055, -0.281, 0.527, 0.246, 0.762, 0.853, 0.891, -0.231, 0.384, 0.373, -0.717, -0.957, -0.338 };
  1621.    double X[] = { 0.884, 0.636, 0.921, 0.282, -0.649, 0.188 };
  1622.    int incX = -1;
  1623.    double x_expected[] = { 0.884, 0.636, 0.653832, 1.112064, -0.168856, 1.225508 };
  1624.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1625.    {
  1626.      int i;
  1627.      for (i = 0; i < 3; i++) {
  1628.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 965) real");
  1629.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 965) imag");
  1630.      };
  1631.    };
  1632.   };
  1633.  
  1634.  
  1635.   {
  1636.    int order = 101;
  1637.    int trans = 113;
  1638.    int uplo = 121;
  1639.    int diag = 131;
  1640.    int N = 3;
  1641.    int K = 1;
  1642.    int lda = 3;
  1643.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1644.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1645.    int incX = -1;
  1646.    double x_expected[] = { -0.59515, 0.077106, -0.27658, -0.637356, 0.407252, -0.308844 };
  1647.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1648.    {
  1649.      int i;
  1650.      for (i = 0; i < 3; i++) {
  1651.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 966) real");
  1652.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 966) imag");
  1653.      };
  1654.    };
  1655.   };
  1656.  
  1657.  
  1658.   {
  1659.    int order = 101;
  1660.    int trans = 113;
  1661.    int uplo = 121;
  1662.    int diag = 132;
  1663.    int N = 3;
  1664.    int K = 1;
  1665.    int lda = 3;
  1666.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1667.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1668.    int incX = -1;
  1669.    double x_expected[] = { -1.46131, 0.537642, 0.624614, 0.762252, 0.326, 0.428 };
  1670.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1671.    {
  1672.      int i;
  1673.      for (i = 0; i < 3; i++) {
  1674.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 967) real");
  1675.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 967) imag");
  1676.      };
  1677.    };
  1678.   };
  1679.  
  1680.  
  1681.   {
  1682.    int order = 101;
  1683.    int trans = 113;
  1684.    int uplo = 122;
  1685.    int diag = 131;
  1686.    int N = 3;
  1687.    int K = 1;
  1688.    int lda = 3;
  1689.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1690.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1691.    int incX = -1;
  1692.    double x_expected[] = { -0.536274, 0.421806, -0.59515, 0.077106, -0.27658, -0.637356 };
  1693.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1694.    {
  1695.      int i;
  1696.      for (i = 0; i < 3; i++) {
  1697.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 968) real");
  1698.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 968) imag");
  1699.      };
  1700.    };
  1701.   };
  1702.  
  1703.  
  1704.   {
  1705.    int order = 101;
  1706.    int trans = 113;
  1707.    int uplo = 122;
  1708.    int diag = 132;
  1709.    int N = 3;
  1710.    int K = 1;
  1711.    int lda = 3;
  1712.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1713.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1714.    int incX = -1;
  1715.    double x_expected[] = { -0.591, -0.084, 0.98216, 0.400464, 0.131806, -0.026608 };
  1716.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1717.    {
  1718.      int i;
  1719.      for (i = 0; i < 3; i++) {
  1720.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 969) real");
  1721.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 969) imag");
  1722.      };
  1723.    };
  1724.   };
  1725.  
  1726.  
  1727.   {
  1728.    int order = 102;
  1729.    int trans = 113;
  1730.    int uplo = 121;
  1731.    int diag = 131;
  1732.    int N = 3;
  1733.    int K = 1;
  1734.    int lda = 3;
  1735.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1736.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1737.    int incX = -1;
  1738.    double x_expected[] = { -1.68293, 0.796222, -0.96062, 0.415172, -0.082386, -0.182748 };
  1739.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1740.    {
  1741.      int i;
  1742.      for (i = 0; i < 3; i++) {
  1743.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 970) real");
  1744.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 970) imag");
  1745.      };
  1746.    };
  1747.   };
  1748.  
  1749.  
  1750.   {
  1751.    int order = 102;
  1752.    int trans = 113;
  1753.    int uplo = 121;
  1754.    int diag = 132;
  1755.    int N = 3;
  1756.    int K = 1;
  1757.    int lda = 3;
  1758.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1759.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1760.    int incX = -1;
  1761.    double x_expected[] = { -1.737656, 0.290416, 0.61669, 0.73853, 0.326, 0.428 };
  1762.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1763.    {
  1764.      int i;
  1765.      for (i = 0; i < 3; i++) {
  1766.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 971) real");
  1767.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 971) imag");
  1768.      };
  1769.    };
  1770.   };
  1771.  
  1772.  
  1773.   {
  1774.    int order = 102;
  1775.    int trans = 113;
  1776.    int uplo = 122;
  1777.    int diag = 131;
  1778.    int N = 3;
  1779.    int K = 1;
  1780.    int lda = 3;
  1781.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1782.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1783.    int incX = -1;
  1784.    double x_expected[] = { 0.27516, -0.544536, -0.10627, -0.988374, 0.229991, -0.711267 };
  1785.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1786.    {
  1787.      int i;
  1788.      for (i = 0; i < 3; i++) {
  1789.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 972) real");
  1790.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 972) imag");
  1791.      };
  1792.    };
  1793.   };
  1794.  
  1795.  
  1796.   {
  1797.    int order = 102;
  1798.    int trans = 113;
  1799.    int uplo = 122;
  1800.    int diag = 132;
  1801.    int N = 3;
  1802.    int K = 1;
  1803.    int lda = 3;
  1804.    double A[] = { 0.002, 0.95, -0.363, 0.084, -0.646, 0.816, -0.407, 0.099, -0.02, -0.906, -0.874, 0.191, -0.328, -0.968, 0.79, 0.826, -0.795, 0.277 };
  1805.    double X[] = { -0.591, -0.084, 0.707, 0.945, 0.326, 0.428 };
  1806.    int incX = -1;
  1807.    double x_expected[] = { -0.591, -0.084, 0.794924, 0.411234, 0.148739, 0.025577 };
  1808.    cblas_ztbmv(order, uplo, trans, diag, N, K, A, lda, X, incX);
  1809.    {
  1810.      int i;
  1811.      for (i = 0; i < 3; i++) {
  1812.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztbmv(case 973) real");
  1813.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztbmv(case 973) imag");
  1814.      };
  1815.    };
  1816.   };
  1817.  
  1818.  
  1819. }
  1820.